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This  paper  investigates  the  electric  power  extractable  from  an  helicopter  conical  nozzle  equipped  with 
thermoelectrical  modules.  The  thermoelectric  nozzle  is  heated  by  the  final  exhaust  gas  from  helicopter 
turbine  and  cooled  by  oil.  A  computer  model  has  been  developed  to  simulate  the  performance  of  the 
thermoelectric  system.  Results  were  obtained  for  various  operating  conditions  showing  that  the  electrical 
power  produced  in  real  operating  conditions  is  significant  but  currently  insufficient  if  we  consider  the 
weight-to-power  ratio.  The  numerical  model  is  also  used  to  optimize  the  electric  power  showing  a  good 
potential  for  the  future. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  increase  in  the  price  of  crude  oil,  accompanied  by  a  bet¬ 
ter  knowledge  of  environmental  problems  associated  with  global 
warming,  resulted  in  an  upsurge  of  scientific  activity  to  identify  and 
develop  environmentally  friendly  sources  of  electrical  power. 

Previous  generations  of  researchers  and  industrialists  have 
worked  on  production  systems  and  the  use  of  different  energy 
sources.  With  the  change  of  economic,  energy  and  environmental 
constraints,  scientists  are  now  trying  to  rationalize  energy  con¬ 
sumption  by  optimizing  processes  and  their  functioning. 

The  aerospace  industry  is  not  an  exception.  A  revolution  is  going 
through  the  aviation  industry,  faced  with  the  obligation  to  change 
its  methods  of  production  and  consumption  to  reduce  energy  costs 
and  reduce  its  environmental  footprint.  Greenhouse  gas  emissions, 
largely  responsible  for  global  warming,  are  particularly  affected, 
while  the  continued  rise  in  oil  prices  highlights  the  need  for  a  bet¬ 
ter  use  of  energy  resources.  In  addition,  the  electrical  equipments 
which  offer  such  benefits  are  now  required  on  board  involving 
additional  power  generation. 
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Thermoelectricity  has  attracted  increasing  attention  as  a  “green” 
and  flexible  source  of  electricity  able  to  meet  a  wide  range  of  power 
requirements  [1].  Thermoelectric  modules  convert  a  portion  of 
thermal  power  into  electric  power  and  only  require  to  function  the 
presence  of  a  temperature  gradient. 

In  aerospace  applications,  such  generators  have  been  used  reli¬ 
ably  for  over  30  years  of  maintenance-free  operation  in  deep  space 
probes  such  as  for  example  the  Voyager  missions  of  NASA  [2].  Mod¬ 
ern  versions  of  these  generators  are  currently  used  in  aeronautics  to 
feed  sensors  for  purposes  such  as  structural  health  monitoring  [3]. 

Thermoelectric  generators  capturing  waste  heat  could  be  used 
differently  in  the  aircraft  industry.  Electricity  can  be  generated 
with  a  higher  power  generation  by  using  exhaust  heat  produced 
by  the  main  engine.  Many  of  the  other  transportation  sectors  and 
in  particular  the  automobile  one  [2,4-7]  share  common  waste  heat 
technologies  and  consequently  will  share  common  thermoelectric 
solutions.  But  aerospace  engine  technology  along  with  mission  crit¬ 
ical  safety  will  require  a  more  careful  and  lengthy  development 
process. 

The  use  of  direct  thermal  to  electric  energy  conversion  offers 
the  following  advantages:  improve  system  efficiency,  no  extra  fuel 
burnt  and  no  moving  parts  to  produce  electricity,  possibility  to 
operate  over  the  entire  aircraft  flight  envelope,  does  not  affect 
engine  operations. 
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Nomenclature 

A  heat  transfer  coefficient  (W/(m2I<)) 

c  specific  heat  capacity  (J/(kg  I<)) 

D  diameter  (m) 

E  total  energy  (J) 

I  courant  (A) 

I<  thermocouples  thermal  conductance  (W/I<) 

L  length  of  the  conical  nozzle  (m) 

H  enthalpy  (J) 

h  specific  enthalpy  (J/kg) 

M  Mach  number 

rri  mass  flow  rate  (kg/s) 

nTE  number  of  thermocouples 

Nu  Nusselt  number 

p  pressure  (Pa) 

Pr  Prandtl  number 

P0  output  power  (W) 

Q.  rate  of  the  heat  flow  (W) 

fleiec  electrical  resistance  (£2) 

Re  Reynolds  number 

Rioad  load  electrical  resistance  (£2) 

Rth  total  thermal  resistance  (K/W) 

S  area(m2) 

T  temperature  (I<) 

u  velocity  (m/s) 

V  volume  (m3) 

Greek  symbols 

a  Seebeck  coefficient  (V/K) 

X  thermal  conductivity  (W/(m  K)) 

p  density  (kg/m3) 

Subscripts 
C  cold  fluid 

H  hot  fluid 

in  inlet 


TEG  into  thermal  energy  systems,  especially  combined  of  heat  and 
power  production  (CHP). 

Although  the  commercial/practical  applications  of  thermo¬ 
electric  generators  still  remain  marginal  in  electric  generation, 
significant  advances  have  been  made  in  synthesizing  new  materials 
and  fabricating  material  structures  with  improved  thermoelectric 
performance  [9].  Efforts  have  focused  primarily  on  improving  the 
material’s  figure-of-merit,  and  hence  the  conversion  efficiency,  by 
reducing  the  lattice  thermal  conductivity  [12].  Recent  advances 
in  materials  and  materials  processing  as  summarized  in  reviews 
[13,14]  have  led  to  higher  non-dimensional  figure  of  merit  values 
and,  thus,  higher  theoretical  conversion  efficiency  of  TE  modules. 

In  addition  to  the  improvement  of  the  thermoelectric  mate¬ 
rial  and  module,  the  analysis  of  thermoelectric  systems  is  equally 
important  in  designing  a  high-performance  of  these  ones  [15-21]. 
The  main  works  consists  in  predicting  the  performance  of  thermo¬ 
electric  generators  with  particular  configuration  of  heat  exchangers 
and  in  studying  the  effect  of  fluid  flow  rates,  fluid  properties 
and  inlet  temperatures  on  the  power  supplied  by  the  system.  For 
example,  Suzuki  and  Tanaka  [18,19]  studied  thermoelectric  power 
generation  with  multi-panels  or  cylindrical  multi-tubes.  The  out¬ 
put  powers  of  the  proposed  system  were  analytically  deduced  from 
heat  transfer  theory,  and  written  by  non-dimensional  functions  to 
reflect  the  characteristics  of  system  design. 

The  aim  of  the  work  is  to  study  the  feasibility  of  generating 
electricity  by  thermoelectric  modules  in  turbomachinery  for  flight 
applications  and  in  particular  for  an  helicopter.  It  consists  in  partic¬ 
ular  to  evaluate  the  electric  power  extractable  from  thermoelectric 
modules  placed  in  the  conical  nozzles  of  helicopter  turbines.  These 
conical  nozzles  are  heated  by  a  hot  fluid  which  is  the  exhaust  gas 
(whose  thermal  power  is  currently  lost)  and  cooled  by  another 
cold  fluid  which  can  be  oil.  The  heat  transfer  through  the  coni¬ 
cal  nozzle  and  the  temperature  variation  along  the  fluid  path  are 
analyzed  mathematically.  Applying  a  heat  transfer  analysis,  their 
temperature  profiles  and  output  electrical  power  can  be  deter¬ 
mined  numerically  using  the  finite  volume  method.  The  physical 
model  enables  us  to  select  the  best  system,  considering  the  ther¬ 
mophysical  properties  of  the  thermoelectric  materials  and  working 
fluids. 


In  the  field  of  aerospace,  it  is  estimated  that  the  solution  could 
reduce  fuel  consumption  by  0.5%,  which,  if  implemented  solely  in 
the  US,  would  save  passenger  and  cargo  airlines  more  than  $12 
million  every  month,  and  reduce  global  carbon  emissions  by  0.03% 
[8]. 

The  production  of  electric  power  on  aircraft  turbine  is  done 
through  the  alternator  taking  a  part  of  the  mechanical  power 
on  the  shaft  of  the  generator.  It  might  be  interesting  in  order 
to  improve  the  performance  of  turbine  power  generation  to  use 
another  method  of  electricity  generation.  The  goal  is  to  reduce  the 
usage  of  mechanical  power  on  the  shaft  generator  and  the  size  of 
the  alternator. 

The  electric  power  generated  by  thermoelectric  modules 
depends  obviously  on  the  nature  of  the  modules  but  also  on  heat 
transfers  on  both  sides  of  these  modules. 

Even  if  thermoelectric  generator  (TEG)  devices  have  a  low 
energy  conversion  efficiency  of  around  2-5%,  they  are  still  strongly 
advantageous  as  compared  to  conventional  energy  technologies, 
not  only  for  their  well-known  advantages  such  as  high  reliability, 
silence  and  low  environmental  impact  but  also  because  of  their 
capability  of  utilizing  amounts  of  waste  heat  as  an  energy  source  in 
a  simple  and  easy  manner.  Total  reliability  of  thermoelectric  gen¬ 
erators  has  been  demonstrated  in  many  applications  [9-11].  For 
example,  Chen  et  al.  [11]  worked  on  the  overall  conversion  effi¬ 
ciency  improvements  and  economic  benefits  of  the  integration  of 


2.  Physical  model 

The  studied  system  is  a  nozzle  with  a  short  length  and  a  conical 
geometry.  Exhaust  gas  circulate  inside  this  nozzle  with  a  low  reduc¬ 
tion  in  pressure  up  to  the  ambient  atmospheric  pressure.  The  Mach 
number  is  about  0.4-0.5  at  the  nozzle  entry  and  is  reduced  to  0.1  at 
the  outlet.  The  idea  is  to  place  thermoelectric  modules  around  the 
conical  nozzle.  The  modules  will  be  kept  pressed  against  the  nozzle 
through  a  second  wall.  This  one  is  cooled  either  directly  by  the  out¬ 
side  air  (convection)  or  by  other  fluids  such  as  air  or  compressed 
oil  (forced  convection).  The  fuel  was  removed  for  security  reasons. 
So  a  coaxial  heat  exchanger  is  used  around  the  nozzle. 

2.1.  Heat  transfer  analysis 

The  physical  model  is  illustrated  schematically  in  Fig.  1.  As 
already  mentioned,  the  thermoelectric  modules  are  sandwiched 
between  the  hot  exhaust  fluid  and  the  cold  fluid  and  they  are  elec¬ 
trically  insulated  by  ceramic  plates  (see  Fig.  2). 

Divergent  conical  nozzles  consist  for  subsonic  flows  in  slightly 
decreasing  the  velocity  and  minimizing  the  total  pressure  losses. 
The  conical  nozzle  performance,  flow  and  heat  transfer  character¬ 
istics  are  strongly  coupled  so  that  the  knowledge  of  flow  conditions 
is  absolutely  necessary  in  the  interpretation  of  performance  param¬ 
eters. 
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Fig.  1.  Schematic  layout  of  the  model. 


The  heat  applied  at  the  hot  junction  is  scattered  and  lost  at  the 
cold  junction,  corresponding  to  the  temperature  gradient  inside  the 
thermoelectric  module.  The  heat  from  the  hot  fluid  is  transfered  to 
the  panel  surface  by  convection,  to  the  cold  surface  by  heat  conduc¬ 
tion  through  the  solid  (metal,  ceramic  and  TE  pellets)  and  finally  to 
the  cold  fluid  again  by  convection.  The  fluids  are  heated  or  cooled 
along  the  flow  path,  and  their  temperature  profile  through  the  path, 
T(x),  changes  as  a  function  of  position  x. 

In  order  to  simplify  the  physical  model  of  the  thermoelectric 
conical  nozzle,  the  following  assumptions  were  made: 


type  of  flow  are  continuity,  momentum  and  energy  equations: 

9(PhS)  9(phUhS)  _n  m 

dt  +  dx  1  j 

d{puUuS)  t  d{puUuUuS)  _  c3ph 

dt  +  dx  “  *  dx  [Z) 

3(/?h£hS)  3(phWh^h5)  _  Qh 

dt  dx  VH  j 

where  pw ,  Uh,  Ph  and  E H,  respectively  represent  the  density,  veloc¬ 
ity,  pressure  and  total  energy  of  the  hot  fluid.  VH  is  the  volume 
occupied  by  the  hot  fluid,  x  is  the  spatial  coordinate  and  t  is  the 
time.  Qh  is  the  heat  flux  at  the  hot  side.  Furthermore,  the  following 
thermodynamic  relations  hold: 


Hh  =  +  2Uh 

1  2 

£h  =  eH  +  -ufi 
PhHh  =  PhEh  +  Ph 


(4) 

(5) 

(6) 


where  Hh  is  the  total  enthalpy,  h\\  the  specific  enthalpy  and  en  the 
internal  energy.  The  system  of  equations  (l)-(3)  is  closed  by  the 
equation  of  state  for  a  perfect  gas: 


Ph  =  Pur  Th 


(7) 


where  r  is  the  gas  constant  and  TH  is  the  temperature  of  the  hot 
fluid. 


-  The  fluid  flows  are  considered  unidirectional. 

-  The  hot  fluid  (air)  flow  is  compressible  and  Newtonian. 

-  The  cold  fluid  (oil)  flow  is  incompressible  and  Newtonian. 

-  All  the  thermoelectric  modules  are  supposed  to  be  conical  and 
are  composed  of  a  single  layer  of  p-n  junctions  as  illustrated  in 
Fig.  2. 

-  The  thermoelectric  pellets  are  homogeneously  aligned  perpen¬ 
dicular  to  the  heat  flow,  combined  slightly  without  an  open  space, 
and  connected  electrically  in  series. 

-  The  electrical  contact  resistance  between  the  p  and  n  couples  is 
assumed  to  be  negligible. 

-  The  material  properties  for  the  TE  pellets  vary  along  the  length  of 
the  conical  nozzle  with  changes  in  temperature  [22]. 

Based  on  these  assumptions  the  energy  balance  equations  for  the 

hot  fluid  and  the  cold  fluid  are: 

2  A  A.  Hot  fluid 

A  one  dimensional  non-viscous  flow  in  a  conical  nozzle  with 

variable  section  S(x)  is  considered.  The  equations  governing  this 


2 A. 2.  Cold  fluid 
Energy  equation: 


,'dTc  i)Tc  \  ,  3 2Tc  Qc 

PcCp'c(  ar+Ucw )  =lchx^  +  v(; 


(8) 


where  pc ,  cp  C,  uc  and  Tc,  respectively  represent  the  density,  specific 
heat  capacity,  velocity,  and  temperature  of  the  cold  fluid.  Xc  is  the 
thermal  conductivity  of  the  fluid,  Vc  the  volume  occupied  by  the 
cold  fluid  and  Qc  is  the  heat  flux  at  the  cold  side. 

The  initial  and  boundary  conditions  were  chosen  as  follows: 


2A.2A.  Hot  fluid.  Since  the  flow  regime  is  subsonic  (0  <  Mach  <  1 ), 
we  define  the  values  for  the  density  and  velocity  at  the  inlet  face 
[23,24].  At  the  outlet  face  we  impose  the  pressure. 


2 A  .2.2.  Cold  fluid. 

Tc(x,0)  =  Tc,in;  Tc(0,0  =  TCjn;  9Tc^’  -  =  o  (9) 

2.2.  Numerical  strategy 


Fig.  2.  Control  volume  of  TEG. 


The  flow  domain  is  subdivided  into  a  finite  number  of  control 
volumes  with  length  Ax;  for  the  hot  fluid  and  Ax,  =  (cos  )  A  x* 
for  the  cold  fluid  (Fig.  2).  All  the  variables  are  stored  in  the  control 
volume  center  (collocated  arrangement). 

For  the  compressible  flows  (Mach  number,  M>0.2)  it  is  well 
known  that  the  efficiency  and  the  accuracy  of  computed  meth¬ 
ods  deteriorate  drastically  when  M  decreases  below  0.2  [23,24]. 
On  the  other  hand,  when  the  Mach  number  remains  uniformly 
small  (below  0.2),  an  accurate  and  useful  approximation  is  to  con¬ 
sider  the  flow  incompressible  [24].  This  observation  has  led  to  the 
development  of  computing  methods  exclusively  suited  to  incom¬ 
pressible  flows.  No  class  of  methods  is  suitable  for  computing 
flows  in  domains  in  which  incompressible  subregions  as  well  as 
compressible  subregions  occur  simultaneously,  or  for  computing 
weakly  compressible  flows.  For  this,  methods  are  required  with 
uniform  accuracy  and  efficiency  whatever  the  Mach  number  is  [23  ]. 
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We  refer  to  such  unified  methods  for  incompressible  and  compress¬ 
ible  flow  computation  as  Mach-uniform  methods.  In  our  work,  to 
resolve  the  Eqs.  (1  )-(3),  we  have  used  a  coupled  pressure  and  tem¬ 
perature  correction  algorithm.  The  description  of  this  algorithm  is 
largely  detailed  and  presented  in  references  [25,26]. 

The  energy  equation  (8)  is  solved  using  a  finite  volume  solution 
method.  Upwind  scheme  is  used  for  the  convection  term  and  a  fully 
implicit  scheme  for  the  time  integration  [27]. 

Because  the  Eqs.  (1  )-(3)  and  (8)  are  coupled,  the  resulting  alge¬ 
braic  equations  have  been  solved  using  an  iteration  procedure,  at 
each  time  step,  until  the  convergence  has  been  achieved. 


Table  1 

Parameters  for  thermoelectric  power  generation  system. 


Variables 

Values  used  in  this  work 

Conical  nozzle 

Length 

0.32  m 

Thickness 

0.002  m 

Inlet  diameter 

0.236  m 

Outlet  diameter 

0.306  m 

Thermal  sources 

Hot  fluid 

Air  (T^n  =  900  K,  uH,in  =  130  m/s) 

Cold  fluid 

Oil  (T/J1  =  358  K,  the  =  0.21  kg/s) 

2.3.  Thermal  resistances 


beck  effect,  the  rate  of  heat  supply  Qh  and  heat  removal  Qc  in  the 
control  volume  (ij)  are  given  by  [29,30]: 


The  thermal  resistances  are  calculated  by  using  classical  corre¬ 
lations.  and  Rc  are  the  thermal  convection  resistances  of  the  hot 
and  the  cold  fluid,  respectively.  Rmetal  and  Ceramic  are  the  thermal 
conduction  resistance  of  the  metal  and  ceramic  materials,  respec¬ 
tively.  ftcont  is  the  thermal  contact  resistance  and  RTE  is  the  thermal 
conduction  resistance  of  the  TEG. 

n  1  ■  n  1  .  n  =  ■  R 

H  AHiSi  ’  c  AcjSj  ’  TE  2ttXte  Axj  ’  metal 

_  In(r2/r1)  +  ln(r6/r5).  _  ln(r3/r2)  +  ln(r5/r4) 

2 7T /-metal  AXj  2 JTAcerarnjc  AXj 

where  r,  i  and  Aq  are  the  radius  of  the  conical  tube  and  the 
heat  transfer  coefficient  between  the  fluid  and  the  conical  tube, 
respectively  (see  Fig.  2).  5/  and  Sj  are  the  wetted  areas  for  hot  and 
cold  fluid,  respectively.  ATE  is  the  average  heat  conductivity  of  the 
module,  defined  by: 


Ate 


ApSp  +  Ansn 
sp  +  sn 


(10) 


where  Ap,  An,  SP  and  Sn  are  the  heat  conductivity  and  the  cross- 
sectional  area  of  p-  and  n-type  elements,  respectively. 

In  order  to  facilitate  the  parametric  design  optimization,  the 
conical  nozzle  model  implements  empirical  correlations  for  hot 
side  AH)I  and  cold  side  Aq  convective  heat  transfers. 

For  the  cold  fluid  flows  through  an  annular  section,  the  heat 
transfer  coefficient  is  determined  using  the  hydraulic  diameter  Dh 
which  is  defined  as: 


(11) 


where  Ac  and  P  are  respectively  the  flow  cross  sectional  area  and  the 
wetted  perimeter.  It  is  this  diameter  that  will  be  used  in  calculating 
parameters  such  as  Re  and  Nu: 


-  For  turbulent  flow  regime,  which  still  occurs  if  Re  >  2300  the  fol¬ 
lowing  correlation  is  adopted  [28]: 

Nu  =  0.023  Re4/5Pr^3  (12) 

-  For  laminar  flow  regime,  the  Nusselt  number  varies  from  3  to  5 
with  variations  in  hydraulic  diameter  and  the  fluid  temperature. 


Qh.i'j  = 

Qc,i,j  =  nxncpnjE 


<xnH.i.s-^y-+K(THls-Tcls) 
«ITQ',s  +  +  K(Tn.i,s  -  Tcls) 


(13) 

(14) 


where  I  is  the  current  flow  through  a  single  thermocouples,  a ,  I< and 
Reiec  are  respectively,  the  Seebeck  coefficient,  thermal  conductance 
and  electrical  resistance  of  a  single  thermocouple.  nTE  is  the  num¬ 
ber  of  thermocouples  in  the  TE  module,  and  nx  the  number  of 
TE  modules  in  a  circumferential  circulation,  and  the  number  of  TE 
modules  in  the  x  direction,  respectively.  TQ,S  and  THtItS  are  respec¬ 
tively,  the  temperatures  of  the  cold  and  hot  junctions.  In  practice, 
it  is  impossible  to  measure  the  temperature  of  both  the  hot  and  the 
cold  junction  (TCj>s,  7h,i,s)  as  the  p-  and  n-legs  are  interconnected 
by  metal  and  are  thermally  in  parallel  between  two  ceramic  plates. 
However,  it  is  feasible  to  calculate  them  using  the  fluid  tempera¬ 
tures  (Tq-,  Th,z): 

,i  ~  QH,i,j  (Rh  +  ^metal  +  ^ceramic  +  ^cont) 


Tc,i,s  —  TC,i  —  Qc,z,j  {Rc  +  ^metal  +  ^ceramic  +  ^cont) 

The  current  flow  through  the  thermocouples  aligned  in  series  for 
a  given  thermoelectric  generator  is  determined  by  the  number  of 
thermocouples  aligned  in  series  and  the  electrical  resistance  of  cir¬ 
cuit  load.  The  total  voltage  is  the  sum  of  the  voltages  generated  in 
each  control  volume,  which  are  wired  in  series  along  x-coordinate 
direction.  Therefore,  the  current  can  be  written  as: 

f  E^i^HXs  -  TcisfnEThn^ 

Kload  + 

where  n  is  the  total  number  of  the  control  volume  and  Rioad  is 
the  load  electrical  resistance.  In  this  work,  the  load  resistance 
is  selected  to  equal  the  effective  internal  resistance  of  the  ther¬ 
moelectric  modules  so  that  the  maximum  power  output  of  a 
thermoelectric  module  could  be  achieved. 

The  output  power  of  the  thermoelectric  generator  is: 

Po  =  V".  (Qh.1  -  Qcj)  (16) 

z — 1 


For  the  hot  fluid  flow  through  the  inner  section  of  the  conical  noz¬ 
zle,  the  above  correlations  are  used  but  the  hydraulic  diameter  is 
substituted  by  the  inner  diameter  of  the  conical  nozzle. 

2.4.  Thermoelectric  modeling 

Since  a  thermoelectric  converter  is  a  heat  engine  and  like  all 
heat  engines  it  obeys  the  laws  of  thermodynamics.  Taking  into 
account  the  heat  conduction,  the  Joulean  heat  loss  and  the  See- 


3.  Results  and  discussion 

In  the  following  numerical  simulation,  the  thermophysical  val¬ 
ues  of  Bi2Te3  semiconductors  will  be  used  for  the  thermocouples. 
The  thermophysical  properties  of  the  metal  and  ceramic  are  deter¬ 
mined  using  the  temperature  of  the  cold  and  hot  fluids.  Table  1 
shows  the  fluid  values  and  the  dimensions  (real  conditions  of  oper¬ 
ation)  for  the  thermoelectric  conical  nozzle  used  in  this  paper.  The 
simulation  was  carried  out  only  for  counter  flow  type  for  the  conical 
nozzle. 
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Figs.  3  and  4  show  the  temperature  and  the  velocity  variation  of 
the  hot  fluid  through  the  conical  nozzle.  For  subsonic  flows,  increas¬ 
ing  the  cross-sectional  area  causes  the  flow  to  decrease  the  velocity 
and  increase  the  temperature.  The  increase  in  the  hot  fluid  temper¬ 
ature  is  a  direct  result  of  the  energy  conservation.  From  knowledge 
of  the  conservation  of  mass  for  subsonic  flows,  the  test  section  can 
be  designed  to  produce  a  desired  velocity  or  Mach  number  since 
the  velocity  is  a  function  of  the  cross-sectional  area. 

Fig.  5  presents  the  evolution  of  the  temperature  of  the  cold  fluid 
through  the  conical  nozzle.  As  can  be  seen  in  this  figure,  the  tem¬ 
perature  variation  tendency  of  the  cold  fluid  flow  along  the  axis 
of  a  thermoelectric  generator  is  similar  to  that  in  an  ordinary  heat 
exchanger.  Nevertheless,  the  variations  in  temperature  are  almost 
linear.  This  special  feature  in  the  conical  nozzle  of  thermoelectric 
power  generation  is  due  to  the  Seebeck  effect:  a  certain  amount 
of  energy  from  high-temperature  fluid  is  taken  away  as  electric¬ 
ity.  Thus  the  heat  flow  from  the  high  temperature  (hot  fluid)  is  not 
equal  to  the  heat  flow  to  the  low  temperature  (cold  fluid). 


x(m) 


Fig.  4.  Velocity  profile  of  the  hot  fluid  along  the  conical  nozzle. 


Fig.  6  shows  the  temperature  variation  on  the  hot  and  cold  sides 
of  the  thermoelectric  modules  along  the  axis  of  the  conical  noz¬ 
zle  for  counter  flow  type.  We  note  that  the  temperature  on  the 
hot  side  of  thermoelectric  modules  decreases  from  the  inlet  to  the 
outlet  of  the  conical  nozzle.  As  the  value  of  the  hot  fluid  velocity 
decreases  along  the  nozzle  (see  Fig.  4),  the  heat  transfer  coeffi¬ 
cients  between  the  hot  fluid  and  the  hot  side  of  the  thermoelectric 
modules  decreases  along  the  conical  nozzle.  Therefore,  all  this  pro¬ 
cess  results  in  a  change  in  the  temperature  profile  on  the  hot  side 
of  thermoelectric  modules.  On  the  cold  side  of  the  thermoelectric 
modules,  the  temperature  increases  from  the  inlet  to  the  outlet 
of  the  conical  nozzle  for  the  counter  flow  type.  Seebeck  effect  can 
generate  the  thermoelectric  power  depending  on  a  temperature 
difference  between  the  hot  and  cold  junctions.  The  Seebeck  elec¬ 
tromotive  force  is  the  sum  of  multiplication  of  the  relative  Seebeck 
coefficient  and  the  temperature  difference  AT  over  all  the  serial 
connections.  The  problem  to  obtain  the  larger  power  is,  therefore, 
how  to  give  the  larger  AT  to  all  the  thermoelectric  modules  existing 
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Fig.  7.  The  comparaison  of  temperature  difference  across  thermoelectric  modules 
for  the  conditions  of  Table  1. 


Fig.  9.  The  variation  of  power  output  with  the  hot  fluid  inlet  temperature  (Oil:  T™  = 
358 1<,  mc  =  0.21  kg/s). 


between  two  hot  and  cold  fluids.  Fig.  7  shows  the  variation  of  the 
temperature  difference  of  the  thermoelectric  modules.  As  shown 
in  Fig.  7,  for  a  given  inlet  temperature  and  velocity  of  the  fluids, 
we  note  an  increase  in  temperature  difference  through  the  coni¬ 
cal  nozzle.  These  results  show  that  depending  of  their  positions; 
each  module  does  not  supply  the  same  electric  power  which  will 
probably  need  the  use  of  electronic  converters  [31  ]. 

In  general,  the  output  power  depends  mainly  on  the  inlet  tem¬ 
perature  and  the  inlet  velocity  of  the  hot  and  cold  fluids  for  a  given 
thermoelectric  generator.  Fig.  8  illustrates  the  variation  of  the  out¬ 
put  power  with  the  inlet  velocity  of  the  hot  fluid.  We  note  that  the 
output  power  increases  with  the  increase  in  the  inlet  velocity.  This 
result  can  be  explained  by  the  fact  that  the  convective  heat  transfer 
coefficient  is  increased  when  increasing  the  velocity  of  the  hot  fluid 
through  the  conical  nozzle.  This  increase  also  brings  an  increase  in 
the  difference  temperature  through  the  modules  and  therefore,  an 


increase  in  the  output  power.  Fig.  9  presents  the  variation  of  the 
output  power  with  the  inlet  temperature  of  the  hot  fluid  for  counter 
flow  type  when  the  inlet  velocity  of  the  hot  fluid  is  1 30  m/s.  We  note 
that  the  output  power  decreases  with  decreasing  the  inlet  temper¬ 
ature  of  the  hot  fluid  for  a  fixed  inlet  velocity.  This  result  is  due  to 
the  fact  that  the  decrease  in  the  inlet  temperature  of  the  hot  fluid 
leads  to  a  decrease  in  the  average  temperature  difference  through 
the  modules,  and  thus  the  output  power  can  be  increased.  Fig.  10 
illustrates  the  effect  of  the  cold  flow  rate  on  the  output  power  for 
two  inlet  temperatures  of  the  cold  fluid.  It  can  be  observed  that 
for  a  given  inlet  temperatures,  the  output  power  increases  with 
increasing  the  cold  flow  rate.  We  note  that  the  increase  in  the  out¬ 
put  power  is  significant  at  the  beginning,  but  the  increase  trend 
becomes  smaller  as  the  flow  rate  increases,  especially  at  lower  inlet 
temperature  of  the  cold  fluid.  It  is  also  found  that  the  output  power 
decreases  with  increasing  the  inlet  temperature  of  the  cold  fluid  for 
a  fixed  inlet  velocity. 
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Fig.  8.  The  variation  of  power  output  with  the  hot  fluid  inlet  velocity  (Oil:  Pf  = 
358  K,  rhc  —  0.21  kg/s). 


Fig.  10.  The  variation  of  output  power  with  the  cold  flow  rate. 
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The  model  permits  to  evaluate  the  output  power  for  different 
inlet  velocities  and  inlet  temperatures  met  in  real  flight  operating 
conditions. 

4.  Conclusion 

Current  needs  for  fuel  efficient,  low  emission,  power  sources  for 
applications  in  sensors  and  subsystems  in  all-electric  aircraft  and 
on-board  spacecraft  has  brought  forth  a  renewed  interest  in  the 
thermoelectric  technology. 

A  numerical  model  for  a  thermoelectric  generator  installed  at 
the  conical  nozzle  level  of  an  helicopter  was  proposed.  This  model  is 
based  on  one-dimensional  differential  equations  representing  con¬ 
servation  equations  for  the  compressible  hot  fluid  (exhaust  gas)  and 
the  incompressible  cold  fluid.  These  equations  are  restructured  and 
linked  to  the  formulations  of  thermoelectric  modules.  The  proposed 
model  gives  much  more  detailed  predictions  for  the  temperature 
variations  through  the  thermoelectric  modules  along  the  conical 
nozzle. 

It  is  found  that  the  output  electrical  power  is  sensitive  to  the 
physical  properties  of  the  hot  fluid  at  the  inlet  face  of  the  conical 
nozzle.  But  the  appropriate  inlet  velocity  should  be  determined  to 
meet  the  optimal  operating  conditions  of  the  helicopter  turbine. 

The  results  of  this  study  show  that  the  electrical  power  produced 
in  real  operating  conditions  is  significant  but  currently  insufficient 
and  in  particular  if  we  consider  the  electrical  power  obtained  per 
unit  of  added  mass  (the  weight-to-power  ratio). 

But  we  may  hope  thermoelectric  device  improvements  to  have  a 
better  efficiency.  For  the  moment,  available  thermoelectric  devices 
on  the  market  have  a  relatively  low  efficiency  and  generally  cannot 
compete  with  alternative  approaches. 

These  powers  can  be  increased  first  by  improving  the  exchange 
between  the  fluids,  further  reducing  the  thermal  resistance  on  the 
cold  side  (increase  in  the  flow  rate,  reducing  the  flow  area  of  the 
cold  fluid). 

It  is  also  necessary  to  improve  the  occupation  of  the  modules 
because  completely  cover  the  nozzle  to  obtain  the  greatest  power 
is  not  the  better  solution. 

Further  work  could  be  done  on  the  simulation  of  a  thermo¬ 
electric  generator  (generator  and  heat  exchangers)  which  can  be 
adapted  to  different  places  to  operate  over  the  entire  helicopter 
flight  envelope. 

Finally,  it  is  important  to  emphasize  that  all  calculations  per¬ 
formed  by  simulation  need  an  experimental  approach  to  fully 
validate  the  results.  We  project  to  carry  out  an  experimental  study. 


In  conclusion,  the  key  challenges  to  push  forward  this  tech¬ 
nology  include:  high  ZT  materials,  high  temperature  modules,  and 
effective  heat  exchangers. 
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